Comparative transcriptomic and metabolomic analyses reveal differences in flavonoid biosynthesis between PCNA and PCA persimmon fruit

The fruit of the persimmon (Diospyros kaki.) has high economic and nutritional value and is rich in flavonoids. Flavonoids are essential secondary metabolisms in plants. The association between persimmon astringency and changes in the proanthocyanidins (a flavonoid subclass) content is well-known. However, information on the relationships between different astringency types and other flavonoid subclasses and biosynthetic genes is more limited. In this study, an initial correlation analysis between total flavonoids and fruit astringency type, and KEGG analysis of metabolites showed that flavonoid-related pathways were linked to differences between mature pollination-constant non-astringent (PCNA) varieties (‘Jiro’ and ‘Yohou’) and pollination-constant astringent (PCA) fruit varieties (‘Zhongshi5’ and ‘Huojing’). Based on these findings, variations in the expression of genes and metabolites associated with flavonoid biosynthesis were investigated between typical PCNA (‘Jiro’) and PCA (‘Huojing’) persimmons during fruit development. The flavonoid concentration in ‘Huojing’ fruit was significantly higher than that of ‘Jiro’ fruit, especially, in levels of proanthocyanin precursor epicatechin and anthocyanin cyanidin derivatives. Combined WGCNA and KEGG analyses showed that genes such as PAL, C4H, CHI, CHS, F3H, F3’5’H, FLS, DFR, ANR, ANS, and UF3GT in the phenylpropanoid and flavonoid biosynthesis pathways may be significant factors impacting the proanthocyanin precursor and anthocyanin contents. Moreover, interactions between the R2R3MYB (evm.TU.contig7272.598) and WD40 (evm.TU.contig3208.5) transcription factors were found to be associated with the above structural genes. These findings provide essential information on flavonoid biosynthesis and its regulation in the persimmon and lay a foundation for further investigation into how astringency types affect flavor components in PCNA and PCA persimmons.


Introduction
Flavonoids are essential secondary metabolites in plants and include more than 10 000 structural variants (Quideau et al., 2011;Mikhailova et al., 2022;Tahmaz and Soylemezoglu, 2022). Flavonoids are divided into six subclasses according to substitutions and B-ring attachments to the basic skeletal structure, namely, flavones, flavonols, anthocyanins, flavanols, flavanones, and isoflavones (Winkel-Shirley, 2001;Sasaki and Nakayama, 2015). These compounds have essential physiological and ecological functions in regulating plant growth and development, flower coloring, fruit flavor, physiological activities, and adaptation to abiotic stress in plants (Havsteen, 2002;Tohge et al., 2018). Flavonoids produced by plants have many health benefits for humans, including antibacterial, antiparasitic, antiinflammatory, anticancer, and anti-aging properties (Dias et al., 2021). Flavonoids have various pharmaceutical activities and often act as antioxidants according to their free radical-scavenging abilities (Kumar and Pandey, 2013). Thus, flavonoids have gained increasing attention and are widely used in the food, cosmetic, and pharmaceutical industries.
Persimmon (Diospyros kaki Thunb.) is a fruit tree that belongs to the family Ebenaceae and has a long history of cultivation (Saleem et al., 2022). As a major fruit variety with a unique flavor, the persimmon has become increasingly popular and has high commercial value in Asian countries . Flavonoids produced by plants have many health benefits for humans and play crucial roles in both the fruit quality and its economic value (Dias et al., 2021;Xie et al., 2022). Persimmons are rich in phytochemicals such as flavonoids, carotenoids, triterpenoids, fatty acids, and vitamin C (Direito et al., 2021), with flavonoids being the main active antioxidants (Sun et al., 2011). The persimmon fruit has a variety of pharmacological poroperties, including antiadipogenic, hypocholesterolemic, antioxidant, anti-inflammatory, and antitumor poroperties, due to its flavonoid components such as hesperidin, naringin, and nobiletin (Direito et al., 2021). Besides, flavonoids such as proanthocyanidins are associated with the astringency and flavor of the persimmon fruits .
Proanthocyanidins (PAs), also known as condensed tannins, are a subclass of flavonoids and consist of oligomers of catechins that are biosynthesized through the flavonoid branch of the phenylpropanoid pathway (Dixon et al., 2005;Zheng et al., 2021). High concentrations of insoluble PAs usually lead to astringency in persimmon fruit. Based on the fruit characteristics, persimmons can be classified into four astringency types, namely, the pollinationvariant non-astringent (PVNA), PCNA, PCA, and pollinationvariant astringent (PVA) types (Yonemori et al., 2000). In China, almost all persimmon cultivars belong to the PCA type, and no PVA and PVNA types are found (Du et al., 2009). The quality and flavor of the persimmon fruit vary greatly, with significant differences between PCNA and PCA persimmons. The fruit astringency type is not only affected the proanthocyanidins content but also the accumulation of total soluble solids, individual sugars, total phenolics, and total flavonoids (Novillo et al., 2016;Yildiz and Kaplankiran, 2018). The influence of the astringency types on variations in the proanthocyanidins concentration is well-known in persimmon fruit (Akagi et al., 2009). However, little is known about the impact of the fruit astringency type on flavonoid metabolic pathway and its associated enzymes, genes, and TFs in persimmon fruit.
This study conducted transcriptomic and quasi-targeted metabolomic analyses to elucidate both gene expression and metabolite accumulated profiles in different stages of PCNA and PCA persimmon fruit. Specifically, dynamic changes in the expression of genes and TFs in flavonoid biosynthesis and the accumulation of a set of flavonoids were analyzed to clarify and compare the secondary metabolism of the persimmon fruit and its complex effects on the astringency and flavor components between PCNA and PCA persimmons. stage (fully developed fruit color without astringency, stage 4). Persimmon fruit from the 'Zhongshi No.5' and 'Youhou' varieties were also harvested at stage 4 and termed PCA1 and PCNA2, respectively. The four development stages (from stage 1 to stage 4) of the 'Jiro' and 'Huojing' fruit were labeled as Jiro_S1, Jiro_S2, Jiro_S3, Jiro_S4 and Huojing_S1, Huojing_S2, Huojing_S3, and Huojing_S4, respectively. Furthermore, fruits of 141 persimmon germplasms were also harvested at the mature stage (stage 4), and these samples were used for total flavonoids content detection. The fresh fruit was frozen immediately in liquid nitrogen and stored at -80°until used for RNA extraction and metabolic analyses.

Extraction and determination of total flavonoids
The total flavonoids of persimmon fruit were extracted and detected according to Han et al. (Han et al., 2021) with a few modifications. Briefly, fruit powder (precise weight 5 g) was extracted with 60% (v/v) ethanol in an ultrasonic bath (30 min). Total flavonoid content was determined by AlCl 3 -(HAc-NaAC) colorimetric method, and rutin with purity = 98% (Solarbio Science & Technology Co., Ltd.) were used as a standard. The absorbance was determined at 420 nm wavelength in a UV spectrophotometer. Total flavonoids content and astringency type of 141 persimmon germplasms are listed in Supplementary Table 1.

Transcriptome data analysis
Total RNA was extracted using the TRIzol Total RNA Isolation Kit (Sangon, Shanghai, China), and a library was established. Bioanalyzer 2100 was used to assess the RNA integrity. The NovaSeq platform (Illumina, San Diego, CA, USA) sequencing generates 150 bp paired-end readings. HISAT2 software  was used to map the filtered reads to the D. kaki reference genome (unpublished). The prediction of new transcripts was performed using StringTie (Shumate et al., 2022). FeatureCounts (Liao et al., 2014) was used to count the read numbers mapped to each gene. The FPKM of each gene was calculated according to the length of the gene and the read count mapped to the gene. DEseq2 (Love et al., 2014) was used to detect the differentially expressed genes (DEGs) between the two groups, with a |log2-fold change| ≥ 1 and padj ≤ 0.05. ClusterProfiler 4.0 (Wu et al., 2021) was used for DEGs in Gene Ontology (GO) (Ashburner et al., 2000) and KEGG (Kanehisa et al., 2004) functional enrichment analyses. Heatmaps and K-means clustering were prepared using the online software Hiplot . The weighted gene coexpression network analysis (WGCNA) were constructed using all genes and were analyzed using WGCNA R package (Langfelder and Horvath, 2008). The networks were visualized using Cytoscape v3.9.1 (Shannon et al., 2003).

Quantitative RT-PCR analysis
The cDNA was synthesized from the high-quality total RNA using TRUE-script First-Strand cDNA Synthesis Kit (Kemix, Beijing, China). Reactions were performed with LightCycler 480 II (Roche), and PCR conditions were 95°C for 3 min, 45 cycles of 95°C for 5 s, and 55-60°C for 30 s. All analyses were conducted with three biological replicates. The relative expression of each sample was calculated by the 2−DDCt method. The persimmon GAPDH gene was used as a reference gene (Du et al., 2019). All gene primers are listed in Supplementary Table 2. 3 Results

Correlation analysis between total flavonoids content and astringency type
The correlation between flavonoids content and fruit astringency were detected in the mature fruit of a natural population with different persimmon cultivars (Supplementary Table 1). The Pearson correlation coefficient (r) between total flavonoids content of persimmon fruit and astringency types of 141 persimmon resources was 0.415**, indicating a significant correlation between total flavonoids content and fruit astringency type (P < 0.01). The above information indicated that total flavonoids were differences between PCNA and PCA fruit.

Comparison of metabolites between mature PCNA and PCA persimmon fruit
To comprehensively define and compare the metabolite profiles between the PCNA and PCA persimmons, we evaluated metabolite compositions using a Quasi-Targeted metabolome. Mature fruits from the 'Jiro', 'Youhou', 'Zhongshi No.5', and 'Huojing' persimmon varieties were labeled as PCNA1, PCNA2, PCA1, and PCA2, respectively. Principal component analysis based on the data for all compounds separated all samples into four distinct groups, with each sample and its replicates forming a separate group, indicating that there were good correlations within group replicates and differences among the different groups ( Figure S1). A total of 889 metabolites were identified in PCNA and PCA fruit, including five metabolite categories, namely, amino acids and their derivatives (170), flavonoids (135), carbohydrates and their derivatives (76), nucleotides and their derivates (65), and organic acid and its derivatives (62) (Supplementary Table 3).
To systematically identify and compare the metabolic pathways between the fruit of the PCNA and PCA genotypes, KEGG enrichment analysis was conducted on the differential metabolites of the four groups (PCNA1 vs. PCA1, PCNA1 vs. PCA2, PCNA2 vs. PCA1, and PCNA2 vs. PCA2). Phenylpropanoid biosynthesis (ko00940), phenylalanine metabolism (ko00360), flavonoid biosynthesis (ko00941), flavone and flavonol biosynthesis (ko00944), and anthocyanin biosynthesis (ko00942) were found to be significantly enriched. The KEGG annotation results suggested that flavonoid metabolism-related pathways were involved in the nutritional value and taste differences between PCNA and PCA fruit ( Figures 1A-D). Given the importance of flavonoids to human health, the subsequent investigation were focused on flavonoid synthesis during PCNA and PCA fruit development.

RNA-Seq of PCNA and PCA persimmon developing fruits
To evaluate flavonoid variations between PCNA and PCA persimmons, the typical PCNA type 'Jiro' and PCA type 'Huojing' were selected for investigation. The fruit was harvested at four stages (S1-S4), namely, the young fruit stage (S1), expansion stage (S2), turning stage (S3), and mature stage (S4). After the removal of lowquality, poly-N, and adaptor sequences, the RNA-seq of the 'Jiro' and 'Huojing' fruit at the four stages yielded 160.80 GB of clean data. The filtered samples contained nearly 6.70 GB of high-quality data with an average Q30 base percentage of 92.42%. Approximately 85.43% of the reads mapped to the reference D. kaki genome, and 4416 novel genes were also identified. The transcriptome sequencing data were confirmed through qRT-PCR. Seven DEGs were randomly selected from the flavonoid metabolism pathway for qRT-PCR verification ( Figure S2). The expression profiles of these genes were consistent with their FPKM values.
The DEGs were compared using DESeq2 software, and selected DEGs were then analyzed. There were 5009, 4023, 6340, and 10 989 DEGs in the Jiro_S1 vs. Huojing_S1, Jiro_S2 vs. Huojing_S2, Jiro_S3 vs. Huojing_S3, and Jiro_S4 vs. Huojing_S4 groups, respectively ( Figure 2A). Venn diagrams showed that 849 genes were differentially expressed in all the comparison groups, suggesting that these DEGs might perform critical functions in the regulation of flavonoids in PCNA and PCA fruit ( Figure 2B). These DEGs were analyzed by the KEGG database to identify their associated pathways. This showed that the flavonoid biosynthesis (ko00941) and phenylpropanoid biosynthesis (ko00940) pathways were signicificantly enriched in the Jiro_S1 vs. Huojing_S1, Jiro_S2 vs. Huojing_S2, and Jiro_S3 vs. Huojing_S3 comparison groups. These results indicated that flavonoid biosynthesis-related pathways may play essential roles in different stages of PCNA and PCA fruit development ( Figures 2C-F).

Differential gene analysis of flavonoid biosynthesis during the development of PCNA and PCA persimmon fruit
Six expression patterns were generated through trend and clustering analyses of the DEGs in the 'Jiro' and 'Huojing' varieties of the four developmental stages, termed Cluster 1-Cluster 6 ( Figure 3A). The genes in these expression profiles were functionally analyzed by KEGG annotation ( Figure 3B). The expression levels of genes in Cluster 3 were higher in the 'Huojing' fruit than in the 'Jiro' fruit. Furthermore, the expression levels of genes in Cluster 3 decreased gradually with fruit maturation. The KEGG pathway analysis showed that genes in Cluster 3 were mainly involved in several primary metabolic processes, such as photosynthesis (ko00196) and starch and sucrose metabolism (ko00500), and secondary metabolisc processes, such as phenylalanine biosynthesis (ko00400) and flavonoid biosynthesis (ko00941). These results indicated that flavonoid metabolism was involved in the development of 'Jiro' and 'Huojing' fruit.
Genes involved in flavonoid biosynthesis were then selected from Cluster III based on the results of the KEGG analysis.  Table 4). The flavonoid-associated DEGs were more enriched in the Jiro_S1 vs. Huojing_S1 group than in the other groups, and most of the DEGs were decreased ( Figure 3C). The analysis of the gene expression levels found that the numbers of DEGs gradually decreased during the developmental process and were specifically highly expressed in Huojing_S1. These results indicated that flavonoid metabolism occurred predominantly during in the early stage (S1) of fruit development and the expression of flavonoid biosynthesis genes in 'Huojing' was significantly higher than in the 'Jiro'.

Construction of a flavonoid coexpression module during PCNA and PCA fruit development
To further identify the specific genes involved in regulating flavonoid metabolism during the development of PCNA and PCA persimmon fruit, 29 057 genes were used in a WGCNA analysis ( Figure S3). To ensure high-scale independence (near 0.9), the b-value was set at 5 ( Figure S3A). The adjacency and topological overlap matrices were then constructed ( Figures S3B, S3C). Based on average hierarchical clustering and dynamic tree clipping, a total of 42 modules were obtained ( Figure S3D). The expression levels of the MEbrown transcripts were found to be significantly higher in S1 compared with S2, and gene expression was significantly higher in Huojing_S1 than in Jiro_S1 ( Figure S3E). The MEblue module contains 3502 genes and KEGG analysis of genes in these modules showed that the genes in the MEblue module were associated with flavonoid biosynthesis during the fruit development, specifically, phenylalanine biosynthesis (ko00400) and flavonoid biosynthesis (ko00941) ( Figure S3F).
The genes with high connectivity in the MEblue module were further investigated as candidate key genes related to flavonoid metabolism. The top 10% of genes in terms of connectivity were selected as potential hub genes. Of these, 12 hubs were identified as potential regulators  Table 5).
Huojing_S1, and 27 DAFs in Jiro_S1 vs. Huojing_S1. Most of the DAFs were significantly higher in Huojing_S1 than in Jiro_S1, with only Procyanidin B2 and Procyanidin B3 were upregulated in at least two stages. The DAFs in the 'Huojing' variety were stronger than in the 'Jiro', which was consistent with the results of the RNAseq analysis. Further analysis revealed that 8, 4, 12, and 10 flavonoid metabolites were stage-specific for S1, S2, S3, and S4, respectively ( Figures 5A, B). In addition, six metabolites were differentially expressed at all developmental stages, including 3,7dimethoxykaempferol-C-glucoside, corilagin, gallic acid, laricitrin, methyl gallate, and phlorizin. These results further confirmed that there were significant differences in the flavonoids between PCNA and PCA persimmons, and that flavonoid biosynthesis pathways play an essential role in PCNA and PCA fruit development. flavonoid contents during PCNA and PCA persimmon fruit development ( Figure 6). Eleven structural genes and eight flavonoids were mapped to the pathway. Two structural genes (PAL and C4H) participated in the upstream phenylpropanoid pathway and nine structural genes (CHS, 2 CHIs, F3'H, F3'5'H, DFR, ANS, ANR, and UF3GT) participated in the flavonoid biosynthetic pathway. PAL2 catalyzes the transformation of phenylalanine to cinnamic acid and as the expression of C4H was significantly higher in Huojing_S1 than in Jiro_S1, there was a upregulation of r-coumaroyl-CoA, with downregulation of cinnamic acid. Subsequently, a series of flavonoid structural genes, CHS, CHIs, F3'5'H, and F3'H showed significantly higher expression levels in Huojing_S1. Thus, some crucial intermediates such as dihydrokaempferol and dihydromyricetin, produced by the enzymes encoded by these structural genes, accumulated highly in Huojing_S1 than in Jiro_S1. All genes related to anthocyanin biosynthesis, such as F3'5'H, DFR, ANS, ANR, and UF3GT, were significantly downregulated in 'Jiro', resulting in lower accumulation of anthocyanins, such as pelargonidin chloride and cyanidin 3-O-glucoside in 'Jiro' than 'Huojing', which might be associated with the differences in flesh color between 'Jiro' and 'Huojing'. However, the levels of the colorless metabolites catechin (gallocatechin) and proanthocyanidin (procyanidin B2 and procyanidin B3) metabolites were significantly higher in Jiro_S4 than Huojing_S4. These findinge demonstrated the main contributions of these eight DAFs to the differences in flavonoid biosynthesis between PCNA ('Jiro') and PCA ('Huojing'), and the critical regulatory roles of 11 genes associated with flavonoid synthesis were hypothesized.

Discussion
As a fruit tree cultivated worldwide, the persimmon has essential ecological, economic, and social value. Presently, persimmon trees are mainly grown for their fruit, which can be eaten fresh or dried. Persimmon fruit varies greatly in terms of specific qualities, with significant variations observed between PCNA and non-PCNA persimmons (Novillo et al., 2016;Yildiz and Kaplankiran, 2018). Volatile compounds also vary significantly between PCNA and PCA fruit, especially in terms of aldehydes (Elhadi, 2017). In addition to soluble tannin, the astringency type also affected the contents of total phenolics, flavonoids, soluble solids, individual sugars, as well as antioxidant capacity  (Novillo et al., 2016;Yildiz and Kaplankiran, 2018). These studies provide a preliminary survey of how the astringency type affects the quality of the persimmon fruit; however, given the importance of flavonoids, the flavonoid composition and potential mechanisms involved in the regulation of flavonoid biosynthesis between the PCNA and PCA fruit varieties still require clarification.
With the rapid development of transcriptome sequencing, many studies have attempted to elucidate the molecular basis of flavonoid biosynthesis via RNA-seq. The key genes involved in flavonoid biosynthesis were identified by stage-specific transcriptomic analysis in the petals of Camellia nitidissima (Liu et al., 2023). Changes in the key genes and flavonoid metabolites were also investigated using metabolomics and transcriptomics in the developing exocarp and embryo of hickory . Metabolomics examines the overall metabolic profile of plant samples through high-throughput detection and data processing (Foito and Stewart, 2018) and can thus provide a reliable method for investigating compounds contributing to the flavor of the persimmon fruit. GC-MS has been used previously to identify volatile components in persimmon fruit (Besada et al., 2013;Elhadi, 2017). Nineteen polyhydroxyphenols were found to be reduced in CO2-treated fruit using untargeted metabolomics analysis, suggesting that persimmon browning might be caused by phenolic compounds . Differences in metabolites between five Japanese persimmons were investigated using NMR (Ryu et al., 2019). In this study, a comprehensive transcriptomic and metabolite analysis was conducted to determine the differences in flavonoid composition between PCNA and PCA persimmons and identify genes related to flavonoid biosynthesis.
Flavonoid biosynthesis is complicated and diverse and requires the substrates derived from the phenylpropanoid pathway . Phenylpropanoid biosynthesis, flavonoid biosynthesis, flavone and flavonol biosynthesis, phenylalanine metabolism, isoflavone biosynthesis, and anthocyanin biosynthesis were observed to be enriched during the development of hickory fruit . In this study, a full-spectrum metabolomic determination of persimmon fruit was performed using liquid chromatography and triple quadrupole mass spectrometry in the MRM mode. This resulted in the identification of a total of 135 flavonoids, greatly broadening our Weighted gene coexpression network analysis. The linkages between TFs and flavonoid-related structural genes in the MEblue module and the sizes of round rectangles were changed according to gene connectivity.
knowledge of flavonoids in ripening persimmon fruit. KEGG analysis of the differential accumulated metabolites in PCNA and PCA persimmons showed that pathways involved in phenylpropanoid biosynthesis, flavonoid biosynthesis, phenylalanine metabolism, anthocyanin biosynthesis, and flavone and flavonol biosynthesis were significantly enriched, indicating not only the high accumulation of flavonoids in persimmons but also that the flavonoids differed significantly between PCNA and PCA persimmons, which may be responsible for the differences in nutritional value and taste between the PCNA and PCA fruit.
In this study, a total of 70 differential accumulated flavonoids were identified between the fruit of the PCNA ('Jiro') and PCA ('Huojing') persimmon varieties at different developmental stages, of which the top three were flavonoids (36), flavones and flavonols (13), and anthocyanins (8). Twenty-eight DAFs were glycosides, which were mainly glycosylated derivatives of quercetin and cyanidin. Quercetin is mostly present in plants as glycosides and has been reported in foods such as onions, apples, broccoli, and tea, and it also has beneficial effects on health (Zheng et al., 2017). Glycoside modifications enhance the water solubility, structural  complexity, and molecular stability of flavonoids (Bai et al., 2022). Thus, flavonoids with glycosides could play essential roles in plant growth, hormonal balance, and the elimination of toxic endogenous and exogenous substances (Wang et al., 2017;Bai et al., 2022); it has, for instance, been shown that delphinidin and its glycosides enhance plant resistance to a wide range of biotic and abiotic stresses (Silva et al., 2017).
In the early stages of flavonoid biosynthesis, the PCNA ('Jiro') and PCA ('Huojing') cultivars differed in the expression of the C4H (evm.TU.contig22.251), PAL (evm.TU.contig9504.51), CHI (evm.TU.contig3165.103 and evm.TU.contig8036.16), and CHS (evm.TU.contig2115.175) genes. L-phenylalanine is converted to cinnamic acid by PAL, the frist enzyme in the flavonoid biosynthetic pathway (Heldt and Piechulla, 2011). In addition to PAL, C4H, CHS, and 4CL play critical roles in the synthesis of crucial secondary metabolites such as lignin, phenolic acids, coumarin, flavonoids, and anthocyanins Xia et al., 2022). Subsequently, r-Coumaric-CoA produces naringenin, which is catalyzed by CHS and CHI (Yuan et al., 2022). CHS is a key initiating enzyme and forms part of a multi-gene family in most s p e c i e s ( N i e s b a c h -K l ö s g e n e t a l . , 1 9 8 7 ) . C H S (evm.TU.contig2115.175) expression was observed to be significantly lower in cultivar 'Jiro' compared with 'Huojing' at t h e e a r l y s t a g e . T w o D E G s t h a t e n c o d e C H I (evm.TU.contig3165.103 and evm.TU.contig8036.16) were also identified between PCNA and PCA persimmons; these two CHI genes were highly expressed in cultivar 'Huojing' at stage 1. The expression levels of these genes might influence flavonoid metabolism in PCNA and PCA persimmons.
LAR and ANS play essential roles in the synthesis of proanthocyanins and anthocyanins; both are downstream genes of the flavonoid biosynthetic pathway and catalyze leucocyanidin into catechins and cyanidins, respectively (Springob et al., 2003). ANR and LAR encode key enzymes involved in the production of 2,3-cis-flavan-3-ols [(-)-epigallocatechin (ECG), (-)-epicatechin (EC), and (-)-epi-gallocatechin 3-O-gallate (EGCG)] and 2,3trans-flavan-3-ols [(+)-gallocatechin (GC) and (+)-catechin (CA)] respectively (Ikegami et al., 2007). Interestingly, there were no differences in the LAR expression level between the two cultivars but ANR (evm.TU.contig4466.754) was expressed at higher levels in cultivar 'Huojing', a finding similar to previous studies showing that the expression level of DkANR was much higher than that of DkLAR during proanthocyanin accumulation (Akagi et al., 2009). This resulted in a lower epicatechin content in the PCNA persimmon and was one of the main reasons for the reduction in the proanthocyanin contents in PCNA types. Cyanidin and pelargonidin usually provide the red pigment in fruit and flowers (Harborne and Williams, 2000). In this study, UF3GT (evm.TU.contig6534.24) catalyzed the formation of cyanidin 3-Oglucoside, and ANS (evm.TU.contig5828.5) also catalyzed naringin to produce pelargonidin chloride; both two genes were expressed at higher levels in 'Huojing' than in 'Jiro', which might result in the accumulation of less red pigmentation in PCNA persimmons.
Genes involved in flavonoid biosynthesis are mainly regulated by the MYB, bHLH, and WD40 TFs and their MBW complex in plants (Hichri et al., 2011), such as rose (Shen et al., 2019) and pears (Premathilake et al., 2020). R2R3MYB TFs are core members of the MBW complex and are involved in the regulation of flavonoid biosynthesis through binding to the promoter regions of structural genes (Yoshida et al., 2015). In pears, PpMYB17 has been shown to positively regulate flavonoid biosynthesis by activating the structural genes PpCHS, PpCHI, PpF3H, and PpFLS in fruit (Premathilake et al., 2020). In persimmon, the combined action of DkMYB2, DkMYC1, and DkMYB4 (MBW) increases the expression levels of the ANR gene involved in the biosynthesis of the proanthocyanin precursor cis-flavan-3-ols (Gil-Muñoz et al., 2020), supporting the above findings on the structural gene ANR and the cis-flavan-3-ols epicatechin content. Besides, DkMYB14 in the Chinese PCNA (C-PCNA) persimmon was found to suppress proanthocyanin biosynthesis and activate acetaldehyde biosynthesis, resulting in the deastringency of the C-PCNA persimmon fruit . MYB82 is involved in trichome development (Liang et al., 2014) and has potential roles in anthocyanin biosynthesis in Arabidopsis (Yang et al., 2013). In this study, an R2R3MYB (evm.TU.contig7272.598) was identified by WGCNA, which was homologous to AtMYB82 and BrMYB82 (Yang et al., 2013), indicating a potential role of MYB82 in anthocyanin biosynthesis regulation; however, the mechanism remains requires further investigation and confirmation.
In conclusion, a comprehensive metabolomic and transcriptomic analysis of PCNA ('Jiro') and PCA ('Huojing') persimmon fruit was conducted. The concentration of flavonoids in 'Huojing' was found to be significantly higher than in 'Jiro' fruit, especially the concentrations of the proanthocyanin precursor 2,3cis-flavan-3-ols epicatechin and anthocyanin cyanidin derivatives. Combined WGCNA and KEGG analyses showed that genes such as PAL, C4H, CHI, CHS, F3H, F3'5'H, FLS, DFR, ANR, ANS, and UF3GT involved in the phenylpropanoid and flavonoid biosynthetic pathways might be the major factors impacting the proanthocyanin precursor flavan-3-ols and the anthocyanin content. Moreover, the R2R3MYB (evm.TU.contig7272.598) and WD40 (evm.TU.contig3208.5) TFs showed significant connections with the above structural genes. This study provides basic information on flavonoid biosynthesis and regulatory network in persimmon fruit and lays a foundation for ongoing investigations on the influence of astringency types on flavor components in PCNA and PCA persimmon.

Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/, https://www.ncbi.nlm.nih.gov/bioproject/PRJNA910302.

Author contributions
YW and YS carried out the experiments and wrote this manuscript. YW and YS participated in the data analysis. HL, ZW, SD, and PS participated in the sampling and investigation; WH and JF conceived and designed the experiments and helped draft and review the manuscript. All authors contributed to the article and approved the submitted version. This article complies with ethical standards.

Funding
This work was supported by the National Natural Science Foundation of China (32071801) and the National Key R&D Program of China (2019YFD1000600 and 2018YFD1000606).

Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

SUPPLEMENTARY FIGURE 1
Principal component analysis score plot of all metabolites in mature PCNA and PCA persimmon fruit.

SUPPLEMENTARY FIGURE 2
Verification of transcriptomic data by qRT-PCR analysis of the expression of seven genes.